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FIELD OF THE INVENTION 

The present invention pertains to the field of communications. More particularly, 
the present invention relates to identifying sources of interference. 



BACKGROUND OF THE INVENTION 

Communication networks are common. Most communication networks 
experience degradation in transmitted signals. This degradation may be from signal loss 
directly, such as smearing of the signal through the medium, loss of signal strength, etc. 
Another source of degradation is noise. Noise may be wideband, narrowband, Gaussian, 
colored, etc. Another source of signal degradation may be from other signals. Often this 
type of degradation or interference is called crosstalk (als o cross-tal k). 

Crosstalk refers to the case signals become superimposed upon each other. The 
signals may be superimposed by electromagnetic (inductive) and/or electrostatic 
(capacitive) coupling in wireline networks. Signals from adjacent transmitters rnay also 
be superimposed over the air in wireless networks. Also, signals from adjacent frequency 
bands or wavelengths may be superimposed in cable and optical networks respectively. 
Crosstalk may come from a variety of physical sources and/or properties, such as bundles 
of twisted pairs that may be capacitively coupled. In bundles of wires, crosstalk may be 
reduced by the use of shielded cables or increasing the distance between the signal 
carrying lines. In wireless and optical networks, crosstalk may be reduced by increasing 
the transmitter and wavelength spacing respectively. Shielded cables are more expensive 
than twisted pair and so this results in increased cost. Increasing the distance between 
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conductors would result in an increased cable bundle size that may present a space 
problem. Similarly, increasing the distance between transmitters or wavelengths in 
wireless and optical networks reduces the system's capacity. Thus, signal crosstalk is a 
problem because it degrades communications. For this reason, the accurate 
5 characterization of the interfering sources may be useful in the analysis, diagnosis and 
ultimately mitigation of the interference. 

SUMMARY OF THE INVENTION 

A method and apparatus for identification of interference sources are disclosed. 
10 Other features of the present invention will be apparent from the accompanying drawings 
and from the detailed description that follows. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention is illustrated by way of example and not limitation in the 
figures of the accompanying drawings, in which like references indicate similar elements 
and in which: 

5 Figure 1 illustrates an exemplary communication system in which the present 

invention may be practiced; 

Figure 2 is a diagram of a DSL communication system in which the present 
invention may be practiced; 

Figure 3 illustrates a bundle of twisted pairs; 
10 Figure 4 illustrates a flowchart overview in which the present invention may be 

practiced 

Figure 5 illustrates a communication channel model in which the present 
invention may be practiced; 

Figure 6 is a flow diagram of one embodiment of an identification process; 
15 Figure 7 illustrates the generation of the /-th disturber from the j-ih service type 

showing the synchronization sequence and the random data; 

Figure 8 illustrates a service type identifier composed of a resampler, a frame 
averager, a matched filter, and a peak detector; 

Figure 9 shows a block diagram of a frequency zoom in algorithm followed by an 
20 FFT analysis; 

Figure 10 illustrates one embodiment of blind baud rate estimation; 

Figure 1 1 illustrates identification using a sequence of known symbols; 
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Figure 12 illustrates one embodiment of a signal flow of a joint co-channel 
identification and symbol detection architecture based on a batch identification algorithm; 

Figure 13 illustrates one embodiment of a batch identification algorithm; 

Figure 14 illustrates one embodiment of a data-aided adaptive algorithm to track 
time-varying co-channels; and 

Figure 15 illustrates where one embodiment of the present invention may be practiced 
in a DSL modem with crosstalk compensation capability. 
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DETAILED DESCRIPTION 

A method and apparatus for identifying interference sources are described. For 
purposes of discussing and illustrating the invention, several examples will be given in 
the context of a wireline communication system, such as DSL. However, one skilled in 
5 the art will recognize and appreciate that interference, for example, crosstalk is a problem 
in wired and wireless communications and that the techniques disclosed are applicable in 
these areas as well. 

In the following description, for purposes of explanation, numerous specific 
details are set forth in order to provide a thorough understanding of the present invention. 

10 It will be evident, however, to one skilled in the art that the present invention may be 
practiced without these specific details. In some instances, well-known structures and 
devices are shown in block diagram form, rather than in detail, in order to avoid 
obscuring the present invention. These embodiments are described in sufficient detail to 
enable those skilled in the art to practice the invention, and it is to be understood that 

15 other embodiments may be utilized and that logical, mechanical, electrical, and other 
changes may be made without departing from the scope of the present invention. 

Some portions of the detailed descriptions that follow are presented in terms of 
algorithms and symbolic representations of operations on data bits within a computer 
memory. These algorithmic descriptions and representations are the means used by those 

20 skilled in the data processing arts to most effectively convey the substance of their work 
to others skilled in the art. An algorithm is here, and generally, conceived to be a self- 
consistent sequence of acts leading to a desired result. The acts are those requiring 
physical manipulations of physical quantities. Usually, though not necessarily, these 
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quantities take the form of electrical or magnetic signals capable of being stored, 
transferred, combined, compared, and otherwise manipulated. It has proven convenient 
at times, principally for reasons of common usage, to refer to these signals as bits, values, 
elements, symbols, characters, terms, numbers, or the like. 

It should be borne in mind, however, that all of these and similar terms are to be 
associated with the appropriate physical quantities and are merely convenient labels 
applied to these quantities. Unless specifically stated otherwise as apparent from the 
following discussion, it is appreciated that throughout the description, discussions 
utilizing terms such as "processing" or "computing" or "calculating" or "determining" or 
"displaying" or the like, refer to the action and processes of a computer system, or similar 
electronic computing device, that manipulates and transforms data represented as 
physical (electronic) quantities within the computer system's registers and memories into 
other data similarly represented as physical quantities within the computer system 
memories or registers or other such information storage, transmission or display devices. 

The present invention can be implemented by an apparatus for performing the 
operations herein. This apparatus may be specially constructed for the required purposes, 
or it may comprise a general -purpose computer, selectively activated or reconfigured by a 
computer program stored in the computer. Such a computer program may be stored in a 
computer readable storage medium, such as, but not limited to, any type of disk including 
floppy disks, optical disks, CD-ROMs, and magnetic-optical disks, read-only memories 
(ROMs), random access memories (RAMs), EPROMs, EEPROMs, magnetic or optical 
cards, or any type of media suitable for storing electronic instructions, and each coupled 
to a computer system bus. 
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The algorithms and displays presented herein are not inherently related to any 
particular computer or other apparatus. Various general purpose systems may be used 
with programs in accordance with the teachings herein, or it may prove convenient to 
construct more specialized apparatus to perform the required method. For example, any 
5 of the methods according to the present invention can be implemented in hard-wired 
circuitry, by programming a general-purpose processor or by any combination of 
hardware and software. One of skill in the art will immediately appreciate that the 
invention can be practiced with computer system configurations other than those 
described below, including hand-held devices, multiprocessor systems, microprocessor- 

10 based or programmable consumer electronics, DSP devices, network PCs, 

minicomputers, mainframe computers, and the like. The invention can also be practiced 
in distributed computing environments where tasks are performed by remote processing 
devices that are linked through a communications network. The required structure for a 
variety of these systems will appear from the description below. 

15 The methods of the invention may be implemented using computer software. If 

written in a programming language conforming to a recognized standard, sequences of 
instructions designed to implement the methods can be compiled for execution on a 
variety of hardware platforms and for interface to a variety of operating systems. In 
addition, the present invention is not described with reference to any particular 

20 programming language. It will be appreciated that a variety of programming languages 
may be used to implement the teachings of the invention as described herein. 
Furthermore, it is common in the art to speak of software, in one form or another (e.g., 
program, procedure, application...), as taking an action or causing a result. Such 
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expressions are merely a shorthand way of saying that execution of the software by a 
computer causes the processor of the computer to perform an action or produce a result. 

It is to be understood that various terms and techniques are used by those 
knowledgeable in the art to describe communications, protocols, applications, 
5 implementations, mechanisms, etc. One such technique is the description of an 

implementation of a technique in terms of an algorithm or mathematical expression. That 
is, while the technique may be, for example, implemented as executing code on a 
computer, the expression of that technique may be more aptly and succinctly conveyed 
and communicated as a formula, algorithm, or mathematical expression. Thus, one 

10 skilled in the art would recognize a block denoting A+B=C as an additive function whose 
implementation in hardware and/or software would take two inputs (A and B) and 
produce a summation output (C). Thus, the use of formula, algorithm, or mathematical 
expression as descriptions is to be understood as having a physical embodiment in at least 
hardware and/or software (such as a computer system in which the techniques of the 

15 present invention may be practiced as well as implemented as an embodiment). 

A machine-readable medium is understood to include any mechanism for storing or 
transmitting information in a form readable by a machine (e.g., a computer). For example, a 
machine-readable medium includes read only memory (ROM); random access memory 
(RAM); magnetic disk storage media; optical storage media; flash memory devices; electrical, 

20 optical, acoustical or other form of propagated signals (e.g., carrier waves, infrared signals, 
digital signals, etc.); etc. 
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Overview of General Communication Network 

The present invention is applicable to a variety of communication systems, for 
example: wireline, wireless, cable, and optical. Figure 1 illustrates an exemplary 
communication system 105 that may benefit from the present invention. The backbone 
5 network 120 is generally accessed by a user through a multitude of access multiplexers 
130 such as: base stations, DSLAMs (DSL Access Mulitplexers), or switchboards. The 
access multiplexers 130 communicate with the network users. The user equipment 140 
exchanges user information, such as user data and management data, with the access 
multiplexer 130 in a downstream and upstream fashion. The upstream data transmission 

10 is initiated at the user equipment 140 such that the user data is transmitted from the user 
equipment 140 to the access multiplexer 130. Conversely, the downstream data is 
transmitted from the access multiplexer 130 to the user equipment 140. User equipment 
140 may consist of various types of receivers that contain modems such as: cable 
modems, DSL modems, and wireless modems. In this network access system the current 

15 invention may be practiced to identify sources of interference in the access channels. 

For illustration purposes and in order not to obscure the present invention, an 
example of a communication system that may implement the present invention, in one 
embodiment, is given in the area of DSL communication systems. As such, the following 
discussion, including Figure 2, is useful to provide a general overview of the present 

20 invention and how the invention interacts with the architecture of the DSL system. 

Overview of DSL Example 

DSL is to be understood to refer to a variety of Digital Subscriber Line (DSL) 
standards that, even now, are evolving. Each DSL standard will be referred to as a DSL ' 
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service type. At the present time, DSL service types include, but are not limited to, 
ADSL, SDSL, HDSL, and VDSL (Asymmetrical, Symmetrical, High speed, and Very 
high speed DSL respectively). 

Figure 2 illustrates a communication system 200, in which the present invention 
5 may be practiced. A central office 202 has a series of DSL modems 204-1 through 204- 
N connected via twisted pairs 206-1 through 206-N as a bundle 208 connected to 
customers DSL 210-1 through 210-N which is connected respectively to customer's 
premise equipment (CPE) 212-1 through 212-N, such as computers. One skilled in the 
art recognizes that twisted pair bundle 208 may experience crosstalk between the twisted 

10 pairs 206-1 through 206-N and depending upon the servic es carried by p airs, data rates, 
an d othe rfactors, such as jro\imit>^ ofthg^ pairs to each other, etc., there may be varying 
a nd different amounts of crosstalk on pairs. 

For example, Figure 3 illustrates a bundle (also called a binder) 308, having 
twisted pairs 306-1 through 306-N. Pair 306-1 may be expected to experience more 

15 crosstalk from a pair 306-2 closer to it than mor e distant 306-L . Likewise, pair 306-2 
located on the perimeter of the bundle 308 may experience different crosstalk than a pair 
306-M more toward the center of the bundle 308. Additionally, if pair 306-1 was the 
only DSL service pair and now pair 306-M is placed into DSL service, there may be new 
crosstalk due to this activation. Also the type of DSL service (i.e. SDSL, etc.) may have 

20 an effect on crosstalk. In general, each DSL service type occupies a band limited 
frequency region. If pairs in proximity to each other are conveying information in 
different frequency bands , then there may be less crosstalk than if pairs are conveying 
information in the same fre quency band . For purposes of discussion, co-channel is used 
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to describe the physical coupling between two interfering pairs. This coupling may be 
represented by a linear dynamic system that will also be called a co-channel. 

Figure 4 illustrates a flowchart overview in which the present invention may be 
practiced. A crosstalk identification device 400, initially acquires signals at 410 that will 
5 be analyzed. At 420 an identification of the crosstalk sources is made and a list of 
models 430 is obtained. At this point, the information may either be stored for later 
analysis or passed onto, for example, another processing step. For example, if the 
purpose of the identification procedure is to enable a crosstalk compensation device, then 
the information may be passed to a compensation design block. It is to be understood 

10 that depending upon shifts, drifts, changes in the communication channel, changes in the 
communications deployed, changes in communications setups, etc., that for optimum 
compensation the steps as detailed above for Figure 4 may be repeated at some interval. 

In order to illustrate the present invention, as mentioned above, the use in an DSL 
system will be described and discussed, however as also mentioned above, one is to 

15 understand that one of ordinary skill in the art will recognize that the techniques 

presented are not limited to DSL and may be used in all manner of communication both 
wired and wireless. 

Description of Received Signal 

20 In order to fully describe the present invention techniques, details relating to the 

signal received at the input of the modem will be described. While one skilled in the art 
will consider this a review, it affords the reader use of a consistent terminology and 
symbol usage for denoting aspects of the invention. The structure of the received signal is 
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depicted in Figure 5 where it is denoted by y(t). In one embodiment of the invention the 
received signal y(t) is sampled by an analog-to-digital converter (ADC) block 540 
producing y(n). We use the notation y(t) to represent continuous time signals and y(n) to 
represent discrete time signals. The discrete time signal y(n) is then passed on to the ID 
5 module for further processing. 

To facilitate the description of the invention, we will focus our explanation on a 
^^^^^^i£iy^P£^tii^^ I" particular, we will concentrate on pulse amplitude modulation 
(PAM) crosstalk disturbers. However, it will be clear to one skilled in the art that the same 

procedure can be applied to other crosstalk sources and is not limited to any particular 

c 

10 modulation technique. Crosstalk sources such as quadrature amplitude modulation (QAM), 
carrierless amplitude and phase modulation (CAP), etc, or any mixture of modulations may 
also be analyzed through the same procedure. 

The received signal y(t) generally consists of a large number of components 
contributed from various sources of signal and interference. Figure 5 describes those 
15 components in more detail. Generally, the received continuous time signal y{t) is 

>'(0 = y„„,>,(0+3'rf,„(0 

where the signal y^Jt) contains the contribution of all the possible disturbers. We will refer to 
this signal as the aggregated disturbance signal. The aggregated disturbance signal can be 
decomposed into two terms: y^Jt) contains the contribution of the PAM signals only, and 
20 v( ') represents the unmodeled noise. 
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Of particular interest is the signal y^^Jt) which contains the disturber signals which we 
wish to characterize. Assume that there are J PAM disturbers that are explicitly modeled in 
the received signal. Then 

y,.nit) = Y^\^^y,{t) (2) 
5 and each individual PAM signal is 

y;(0= X^,W^(^-A:r.) (3) 
where Sj{ k) represents the transmitted PAM sequence of the 7-th disturber through an overall 
co-channel impulse response hj{t) and with symbol period Ty. Finally, the received signal 
sampled at a sampling rate T is 

10 ypoM = ypa.it)[._,,^ -^yjin) (4) 

where 

y>(«) = S^; W^ (^^. -^^y) (5) 

it 

The noise signal has little structure and is the simplest of the three components to 
describe. The sampled noise signal 

15 v(n) = v(0^ (6) 

may be modeled as additive Gaussian noise the color of which is characterized by the power 
spectral density of the signal v(t). In other cases, the noise term can model other interfering 
signals that will not be actively characterized like impulsive noise, AM radio interference etc. 
Last, but not least, the main signal may be present in the received signal. 

20 If the service type on the main line is the same as the service type on the disturber lines, 

I 
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then the main signal will have an identical description with the one given above for each 
disturber signal. If the service type on the main line is for example an ADSL service, 
then the main signal will employ DMT nriodulation and its description will be different 
(for details on the description of a DMT main signal see co-pending patent application 

Serial No. titled "Method and Apparatus for Mitigation of Disturbers in 

Communication Systems" assigned to the assignee herein and filed on even date 
herewith. In several cases the receiving modem may be able to force the modem 
transmitting on the same line to silence through the use of an appropriate command. In 
that situation there is no main signal component in the received signal and y^^.^ («) = 0 . 
In any case, the main signal does not play an important role in the disturber 
characterization process and its exact description is not required in the current context to 
understand the present invention. In fact, the main signal is removed from the received 
signal before the identification proceeds as described next. 
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Signal Characterization Procedure 

Now that the received signal has been described, an overview of the steps that 
occur during identification time will be discussed. There are four main steps: 

1. detection of service types present, 

2. baud rate estimation, 

3. setup of co-channel identification, and 

4. initial co-channel identification. 

An overview of each process will be given, with details to follow. 
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Detection _of service types present is a technique that determines the frequency 
regions with significant disturber energy. Since there are a large number of possible baud 
rates that may impair the main line, it is not realistically feasible at this time to try each 
single rate in order to determine if a service type is present or not. Therefore, an initial 
5 coarse selection of the possible frequency regions containing disturber energy accelerates 
the entire identification process. The outcome of this process is a collection of data rates 
which contribute disturber energy to the received signal. 

Each data rate determined from the above process represents a possible disturber 
service type present in the transmission. However, several disturbers may correspond to 
10 the_saine_se^^ Once the possible disturber data rates and/or 

service types are determined, the accurate^baud rate and co-chan nel est imation steps are 
repeated for each identified service type. 

Note that due to oscillator differences, the actual timing signal used by the 
disturber generation may not be synchronized with the main channel timing signal. For 
15 example, crystal oscillators are known to differ from the nominal frequency by as much 
as 100 parts per million. However, in order to obtain acoi rate co-channel estimation, we 
estimate the difference between the ^ctual^PAI^^ the nomina[^ne This is 

done during the baud rate estimation step. 

Figure 6 is a flow diagram of the overall identification process. The first step in 
20 the process for the present invention in one embodiment in a DSL modem would be the 
collectioi^^ Note that with specific reference to an 

ADSL modem, the identification operations may be performed during Medley, after time 
equalization (TEQ) and frequency equalization (FEQ) training. Thus, in order to obtain 
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the aggregate disturbance signal, one would need to remove the main signal. For an 
example of particular details of the main signal removal procedure see co-pending patent 

application Serial No. titled "Method and Apparatus for Mitigation of 

Disturbers in Communication Systems" assigned to the assignee herein and filed on even 
5 date herewith. It may be possible that in other uses of the invention, the signal from the 
main channel is not present during identification time. This may happen for example if 
identification is performed before the main channel transmitter is powered on or is 
otherwise allowed to transmit, or is instructed to not transmit. Then, the received signal 
is simply the^ggregated disturbance si gna|. It is clear that in this situation the main 

10 signal removal step is not required. 

The next step during identification is the detection of the service types present 
604 in the signal. Next is a sequence of three major steps that may be repeated for each 
service type present 610. The first step in the three major steps is that of a baud rate 
estimation 606, followed by the second step, a setup of the co-channel identification 

15 procedure 607, and the third step is an initial co-channel identification 608 using symbols 
embedded in the signal that are known a priori. The result after step 608 is an initial 
model of the co-channel. If more service types remain unprocessed, then for each service 
type present 610 the steps 606, 607, and 608 are repeated. When all service types present 
have been processed the result is a list of models 612. This list of models 612 may be 

20 used to create, construct, modify, and/or design a compensation system. Alternatively, 
the list of models niay be used to analyze the crosstalk disturbance of a particular 
communication channel. 
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The list of models 612 can be further refined during a final co-channel estimation 
614. Finally, when dealing with long observation periods, time varying co-channels, etc., 
a parameter adaptation procedure 616 may be advantageous to implement. 

Next, we will describe parts of the identification process in more detail. 

5 

Service Type Identification and Baud Rate Estimation 

Existing baud rate estimation techniques for single disturbers use nonlinearities to 
obtain a periodic signal with a period that is the desired baud rate, and then use a Phase Lock 
Loop (PLL) to track small differences. However these techniques require good SNR levels in 

10 order to detect small phase errors, and they cannot be applied when several signals with 
similar baud rates and energy levels are present at the same time. To circumvent some of 
these problems, the present invention exploits the cyclostationary properties of the disturbers 
and performs a search in the frequency domairL The resulting technique is accurate and may 
be implemented in an efficient form. 

15 To start the description of the procedure, let us rewrite Equation (1) that describes 

the aggregated disturbance as follows: 

y^is^ («) = S Z S 'ji (''^'^ji - ^7^.7 ) + Kn) (7) 

j I k 

In this equation, the index ; goes through the set of service types, / indexes among 
20 all the disturbers from the j-th type, s/k) is the sequence of symbols sent by the l-ih 

disturber of type/ Similarly, is the baud rate, and is the co-channel for the l-ih 
disturber of type y. 

Note that the actual baud rate T.^ may have an offset with respect to the ;-th 
nominal frequency. This offset is determined by the characteristics of the local oscillator 
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in the disturber transmitter. The local oscillator at the disturber transmitter determines 
the actual baud rate of a particular disturber, as well as its timing signal. In general, the 
local oscillator has a constant unknown offset with respect to its nominal frequency that 
can cause maximum frequency errors of 100 parts-per-million. The maximum allowable 
5 frequency offset for a particular disturber type is specified in the corresponding service 
type standard. If the observation time is short enough, it is possible to neglect 
instantaneous phase errors of the timing signal due to frequency drift with time and other 
random effects. The use of a short segment is advantageous from an implementation 
perspective, and under these conditions we may assume that the only source of phase 

10 error is a constant frequency offset with respect to the nominal frequency. The issue of 
timing signal tracking for longer periods is considered in a later section below. 

Notice that the received disturber signal is a mixture of transmitted signals of 
diffe r ba^ idjates. In several applications, the nominal frequencies of the disturbers 
may be unknown. Even when several disturbers of the same nominal frequency are 

15 pjisent, the actual individual baud.rates may be different due to the differences among 
the local oscillators in the disturber transmitters. 

It is also important to observe that the co-channels h.{.) may have comparable 
energy levels. Therefore, some of the individual disturbers in the received mixture may 
have similar levels of total energy. This implies that in general, the signal to noise ratio 

20 (SNR) of any given disturber computed as the ratio between total signal energy for the 
said disturber and total noise and interference energy may be very poor and traditional 
baud rate estimation techniques may fail in this situation. An alternative approach is to 
perform a precise search in the frequency domain using a Fast Fourier Transform (FFT). 
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To assist the reader in understanding, we will first analyze the case of baud rate 
estimation when a single disturber is present. Let T be its corresponding baud rate. The 
technique described in this section is based on the cyclostationary properties of the signal 
in Equation (7). For that, we estimate the correlation of the signal y^,Jn), Let r (rt,x) be 
the time varying autocorrelation of y ^Jn) as follows 



It is possible to demonstrate that r (n,0) is a periodic signal with period T. Then, 
the baud rate 7 can be recovered as the period of r{n,0). Let us denote by r(n) an 
instantaneous estimate of r(n,0), i.e., 



Then it can be shown that the Fourier transform of r(n) converges to the Fourier 
transform of r(n,0) as n goes to infinity. Hence, r(n) can be used to estimate the period T, 

In the case where multiple disturbers are present, if the disturbers are independent 
one from the others, and the co-channels are different, r{n) will contain the sum of 
processes with periodic components. Therefore a careful search for the periodic 
components of r{n) will yield the desired answer. One possible technique to perform this 
search is to use a fast Fourier transform (FFT). In order to reduce the complexity of the 
overall technique without reducing the accuracy of the estimation, we will select 
candidate frequency regions to perform the searches. It is possible to improve the 
resolution of an FFT along a certain frequency region by "zooming in" the desired 
frequency region. For example, if the desired frequency region is centered oif^ and has a 
bandwidth Wthen it is possible to modulate r{n) by the nominal frequency/^, and the 
resulting signal is r^^(n) 



r(n,T) = £[v (n)v (n + T)] 



(8) 



(9) 



(10) 



00I34O.PO77 



# 



After removing the high frequency components from r ,(n), the resulting signal is a 
baseband signal 



where /7^p(n) is a low pass filter with cutoff frequency/y2. 

It is possible to reduce the bandwidth of r^(n) to be equal to Why using a cascade 
of lowpass and decimator filters. For example, if we let L be the decimating factor and 
h^pj the lowpass filter corresponding to the first decimating filter then, the output after the 
first decimation is as follows 



Notice that the bandwidth of r^^^iin) has been reduced by a factor L. By applying a 
cascade of low pass and decimator filters, it is possible to reduce the bandwidth of the 
signal r^(n) to W, the bandwidth of the desired frequency region. Then a simple FFT 
analysis allows us to obtain all the harmonic components of the signal in the frequency 
range + W], It is clear that this frequency range corresponds to the frequency range 

Figure 9 shows a block diagram of the frequency zoom in algorithm 910 followed 
by a FFT analysis 920. 

From the discussion above, the election of the frequency search regions can be 
seen to be an area of practical concern. The nominal frequency is a characteristic of the 
disturber type and it is specified by the service type standard (SDSL, ISDN, etc.) that 
defines each particular disturber type. On the other hand, the bandwidth W is determined 
by the accuracy of the local oscillator, also specified in the applicable standard. If we let 



(11) 



(12) 



k=nL I 
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N be the number of possible disturber types present in the mixture of Equation (1), and 
assume that for each disturber type the standard specifies n. possible nominal frequencies, 
then the set F defined as 

5 is the set of all possible nominal frequencies. When this set is a reduced set of 

frequencies, then it is possible to specify a reduced set of intervals [f^j - W.,f^j + W^,. J. 

However, in certain applications of the present invention, the set F may be very 
large or even unknown. In these cases, an a priori specification of the search regions is 
unfeasible. Nonetheless, it is always possible to perform a cjoarsej^it^^ to 

10 determine the main frequency regions (as illustrated at step 604 in Figure 6) that contain 
significant energy using the frequency zoom in algorithm described above. Another 
alternative is to divide the total bandwidth of the signal r{t) in regions, each one with 
bandwidth IN and then perform a frequency zoom in and an FFT analysis in each 
region. Those frequency regions that exhibit some periodic energy may be further 

15 refined. This procedure may be iterated several times until the desired accuracy in the 
frequency estimation is obtained. This approach is denoted as blind baud rate estimation 
and is further illustrated in Figure 10. 

In order to characterize a particular disturber it is important not only to determine 
it s baud rate b ut also its alphabet, frame length, etc. This particular characterization is 

20 known as service type identification. Once the service type of each disturber has been 
determined, it is possible to use the known features of each standard to perform co- 
channel identifi cation. For example, most DSL standards specify a sequence of 
sync hronization s ynibols^hat are sent per iodically by the centr al office. The time elapsed 
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betweeo_tw o consecutive sync hroni zation sequences isknown as a frame of data. The 
frame length and the synchronization sequence used by each particular standard may be 
known a priori. As it will be explained later, the sequence of known symbols is used as a 
training sequence during the co-channel identification step. It is clear that any sequence 
of symbols known a priori can be used in this procedure. In particular we will describe 
the procedure using the synchronization symbols embedded in the disturber signal. 

In the discussion of the present invention, so far, we have described use of the 
coarse frgguency^es t^ step as a means to determine the service type of the 
disturbers that are present. It will be clear to one skilled in the art that an alternative 
approach is to use known symbols to pe rf orm service type identification. Please note that 
either of these approaches, and others, may complement or supplement one another. An 
advantage of doing s ervice type identification by u singbaud rate estimation is that it ~ 
requires using only a small segment of data, thereby eliminating distortion of the data 
introduced by phase errors of the timing signals 

Initial Co-channel Identification 

Referring to Figure 6, we will describe the setup of the co-channel identification 607 
procedure. To illustrate the procedure, let us focus again on Equation (7). The sequence s/k) 
is divided in two subsequences:' 

s,ik) = s](k)^s'j^(k) (13) 
The first subsequence s''ji(k) corresponds to the random data; the second subsequence, s)(k) 
corresponds to the known symbols for service type Both sequences are orthogonal, as 
shown in Figure 7, which illustrates the generation of the /-th disturber from the y-th service 
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type. In this example, the synchronization symbols are known. Figure 7 shows the 
s>Tichro^ and the random data. The sequence of known symbols is a periodic 

sequence and its period is the frame length corresponding to the particular service type. 

To illustrate the use of known symbols in co-channel identification, let us assume that 



5 the disturbance signal is re-sampled at a fraction l/P. of the baud rate so that 
y,is. («) = E E Z (^)^7 - I'Pj ) + v{n) (14) 

j I k 

For each service type, we may implement the system shown in Figure 8 that is composed of a 
resampler 801, a f rame averager 802 , a match filter 804, and a peak detector. Let nf. be the 
frame length corresponding to the 7-th service. If we assume that NmJ. samples have already 
10 been collected. Then, the output ofthe frame averager is a sequence ^f^g^n/ that is^ 
obtained as follows: 



3'r(«) = Vt (n + l.nf^) n = 0,-,«/. (15) 

/=0 



In order to determine the presence or absence of the known sequence of symbols, the design 
15 of a matched filter uses the sequence of known symbols, s/0).. . s.(M.-l ), convolved with the 
pulse-shaping filter of the 7-th PAM disturber p. (n) as an approximation to the actual co- 
channel. 

Fj in) = (/) * pj = 5; s] {k)p. ((-« - kP. )T. ) (16) 

k 

Then, 

20 4, ^ (17) 

/ k 



24 



00I340.P077 




When;-th type is present in the mixture of disturbers ( y^.^^ {nT^ ) ), the output of the 7-th 
matched filter has a peak. Peak detection is done using an appropriately selected 
threshold. The value of n corresponding to the peak matches to the position of the center 
of the sequence of known symbols in the averaged fr ame of data. The pe ak detect ion 
5 module generates two important outputs. The first one js th e position of t he 

s ynchronization sequence within the frame of data, which is obtained by observing the 
index n at which a peak is detected. The second output is the n umber of disturb ers of the 
same^tyge^that are present at the same time. This output is obtained by counting the 
^numb^^ 

10 Notice that as was mentioned before, it is possible to use the system described in 

Figure 8 to perform service^p e^dentification. Suppose that the service types that are 
present in y^Jn) are unknown. Then, one can implement a bank of systems such as the 
one shown in Figure 8. Each system will be designed according to a possible service 
type. When a particular service type is present, the output of the matched filter 804 

1 5 displays ^ peak that canjejetect^^ peak detector 806. To decrease the risk of 
false alarm, a hypothesis test may be run with the outputs of the peak detectors. 

Once the ^o-channel identification set up is completed, the averaged frame 
described in Equation (15), the total number of disturbers, and the position of the 
synchronization sequences are passed to the initial co-channel identification step (block^ 

20 608 in Figure 6). 

To illustrate the identification procedure, we decompose the averaged data ( y/"' 
(k) ) into three terms: 
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y'r (") ={y:^> («))+ 'y^s^knf) ^'''^^ 

= )',w,. («) + vv/zi) 




(18) 



where 



y,., («) = in) * hj in) + 2 ^„ in) * h,„ in) (19) 

m 

The first term in Equation (18), y.^^ (n) , corresponds to the contribution of the sequence 

5 of known symbols s.(0). . . sfM.-l) to the y-th disturber. The second term in Equation 
(18), represents the contribution of the random data of the same type s .(k), as well as the 
contribution of the disturbers of different types. Finally, the last term in Equation (18) is 
the noise term. For the purpose of an identification technique, the noise term w.{k) may 

considered to be composed of the contribution of the random data and the additive 

10 Gaussian noise. 

We will use the sequence of known symbols s.(0). . , s](M.-l ) as the input to the 
system. To decre as^ e influence of t hejgndom disturber symbols , we estimate the 
position in the frame of data of the starting time for the sequence of known symbols. Let 
K. be this position, and L be the length of the 7-th co-channel. Then, the identification 

15 problem is to minimize the following cost function among a selected family of models n. 

^=argmin,,n ^ ()'r(^)-^;(^)*/'(^))' (20) 

k=Ki 



The sequence of known symbols is in general very short. For example, in the case 
of HDSL services, only 7 symbols out of 2351 are used for synchronization purposes. 
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Therefore, one cannot assume to have perfect input excitation from this class of signals. 
This also implies that the model order of the models in FI cannot be chosen arbitrarily 
large. 

T o redu cejliehigh ^ noise, both inputs and outputs may be filtered using 

a lowpass as shown in Figure 1 1 illustrating identification using a sequence of known 
symbols. Notice that since both input and output are lowpass filtered, only the noise 
system is modified. However, if only the output data had been filtered, then the 
input/output model, or h(k), would also include the lowpass filter, which may have an 
undesirable effect. 

Several disturbers of the same type may be treated jointly when the starting 
position of their synchronization sequences are close together. In this case, the 
summation as in hj above is extended as follows 

i^,,,.",Vl=argmin ^ br(^)-^; W*^ W) £ {yT(k)'s]ik)^h, 

(21) 

An effective solution for this identification problem is to consider O as a family of 
multiple-input-single-output (MISO) models. Then, standard MISO system identification 
techniques can be applied to this equation ( l^i.^ , • • , } )• 

So far, we have assumed that in order to synchronize the samples in the averaged 
frame, we resample the data using the baud rate estimation. However, when more than 
one disturber of the same type is present this is not possible. Since each disturber 
transmitter may use different oscillators, the actual baud rates of disturbers of the same 
type may differ by small amounts. When the synchronization sequences of these 
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disturbers are fired closely together in the data frame, the contribution of the different 
disturbers cannot be separated, and a MISO system needs to be identified. If the data is 
not resampled at the exact baud rate of a certain disturber, its pulse will be smeared in the 
average frames. Since we know the difference between the resampling rate and the 
5 actual baud rate, it is possible to account for this effect in the identification process. 

For the understanding of the reader, we will present the technique for co-channel 
identification for the case that a single disturber is present in the mixture of disturbers 
ydi.J^)- This result will be extended later to the case of multiple disturbers present in 
y^J^)- L et rbe the baud ra te an d fef^J the co-chaTmel i mpu^ response corresponding to 
10 the single disturber present in y^.^^ (n). In this case. Equation (7) can be re- written as 
follows, 

yi^T^ ) - E ^ik)h{nT^ - kT) + v{nT^ ) (22) 

k 

Using linear interpolation, we express h(nT-kT) as follows: 

h(nT^ - kT) = KOi - k)T)il + n^)- h({n -k- l)T)n ^ (23) 
1 5 In general, using a 21 order interpolation, h(nT-lT) has the following expression: 



h(nT,-kT) = [q^,in,[) ■■■ q,,in,2li 



= q,r{nMin-k)T) 



h{in-k-l)T) 



h{{n-k)T) 



h(in-k + l)T) 



(24) 
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where AT = T-T. Notice that in Equation (22), s(k) is a scalar. Moreover, the vector 
q^r('0 introduced in Equation (24) is independent oik. Thus, q^j.(")can be factored out 
of the convolution summation in Equation (22) as follows 
yinT^ ) = q^, {n)^ s{kM{n - k)T) + v(nr, ) 

k 

= q^T- ('OS - nMkT) + v(nT^ ) 



(25) 



10 



For simplicity, we will develop the procedure for h(.) being an fini te imp ulse response 
(FIR) channel. Howevef, it is straightforward to extend the results from these notes to an 
infinite impulse response (IIR) model. Let L be the length of the co-channel, and H a 
21+L vector constructed from the impulse response h(kT) as follows 



H = 



HO) 
KiL-l)T) 



Then, Equation (25) can be rewritten as follows 

3'(«7:j=q^r(n).^[0 sik-n)l, 0}l + v{nT^) 

n 

= q^r(«)s(«)H + v(n7,) 



(26) 



(27) 



15 



where I, refers to the /-by-/ identity matrix. Now suppose that the data frame is nf 
symbols long. Moreover, suppose that //.«/symbols have been collected. Then, using 
Equation (27), we compute the averaged frame of data as follows 
inT^ ) = i^yin + l.nf) + v((n + l.nf)T^ ) 

= ¥ E^ir (" + lnf)sin + l.nf)H + v((n + l.nf)T^ ) (28) 



^q^r(" + ^"/)s('0 
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The matrix s{n) introduced in Equation (28) contains the known symbols and the random 
data, i.e., 

s(n) = s'{n) + s''{n) (29) 
5 In Equation (29), s' (n) is formed from the sequence of known symbols, and s' («) is 

obtained from the random data. Notice that the sequence s(n) is zero before the first 

known symbol has been sent, and after the last known symbol has been sent. Therefore, 

the structure of s' («) depends on the location of the known sequence within the data 

frame. 

10 We can separate the two components ofs(n) to emphasize the averaging action. 

Lumping together the noise term and the contribution of the random data in a single term 
£, we rewrite Equation (28) as follows: 



s^(n) H + £ 

(30) 



y"-(n7;) = ^|Xq.r(« + ^.«/) 
= Q(n)H + £ 

15 Equation (30) can now be used too btain an FIR model for hLl 

In order to obtain an IIR model, we need to include the contribution of past output 
values in Equation (30). Notice that the interpolation process from Equation (24) can be 
associated to the moving average (MA) portion of an IIR model. Thus, a similar process can be 
applied to the autoregressive (AR) portion of the IIR model. In general, an m-th order IIR 

20 model is expressed as follows 

y{k) -a,y{k-l) y(k -m) = b^s{k) + 6,5(/: - 1) + • - + b^j{k -m) (31) 

In Equation (31), s(k) is the input and y(k) is the output of the system. Let A be the matrix 
formed as in Equation (26) using the coefficients Similarly, we can form the 

matrix B using the MA coefficients Finally, we denote by y'""{k) the vector 
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formed as in Equation (27) using past output values. Then, Equation (30) can be re- 



written as follows: 



(32) 



5 In this case, Equation (32) is the one used to perform co-channel identification. 

We will now discuss the case of multiple disturbers. If several disturbers fire their 
synchronization sequences close together, then we use MISO techniques to obtain the estimates 
of the co-channel. 

When multiple disturbers are present, it is straightforward to pose Equation (30) or (32) 
10 as a multiple input single output system (MISO) 



Notice that in this case the matrix Q/n) is different for each disturber because it 
depends on the location of each synchronization sequence. Co-channel estimates are obtained 
by using standard MISO identification techniques to Equation (33). 

15 The disturber signal y^Jn) in Equation (1) has been described so far as a mixture of 

disturber signals plus additive color noise. The p.urpx)s e of the co ^annel identification 
procedur e^is to describe the stru ctur e of the disturber signah So far we have described how to 
describe the mixture of disturbers. The remaining component of the disturber structure is the 
residual noise temvi^ in Equation (I). To complete the description of the disturber structure, 

20 we obtain a description of the random signal v(t). We will consider that v(t) is a zero-mean 

Gaussian random process. The power spectral density of this signal can be computed using the 



j 



(33) 
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prediction error obtained from the co-channel models previously identified. An example of 



such computation is described in U.S. Patent Application Serial No. 09/523,065, filed March 
10, 2000, entitled Method for Automated System Identification, to C. Galarza, D. Hernandez, 
and M. Erickson, and assigned to Voyan Technology Corporation of Santa Clara, California. 
5 In some applications, it may be necessary to compute the noise model when only one disturber 
co-channel is considered at a time. The noise models obtained with this procedure are 
important for successive disturber cancellation. For an explanation, see co-pending patent 
application Serial No. titled "Method and Apparatus for Mitigation of 



Disturbers in Communication Systems" assigned to the assignee herein and filed on even date 
10 herewith. 

Also in some applications, it is useful to compute uncertainty bounds for the co- 
channel models obtained previously. These bounds are a measure of the identification 
error due to incorrect model structure, finite number of data points, and measurement 
noise, among other factors. It will be appreciated that modeling errors is well known by a 
15 person of ordinary skills in the art and thus further details of this topic will not be 
described. One embodiment of the computation of the uncertainty bounds has been 
described in detail in U.S. Patent Application Serial No. 09/345,640, filed June 30, 1999, 
titled "Model Error Bounds for Identification of Stochastic Models for Control Design", 
assigned to Voyan Technology Corporation of Santa Clara, California. 



Final Co-channel Estimation 

After some initial estimates for the co-channel impulse responses are obtained via 
the method described above, further estimation accuracy may be achieved by employing 




20 
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a data-aided identification procedure described next. This is performed in block number 
614 in Figure 6, which is implemented as an approximation of a joint co-channel 
identification and symbol detection technique. 

One alternative approach is to utilize the batch algorithm approach that will be 
5 described next. To simplify the explanation, the algorithm is described for one type of 
disturber modulation only, namely pulse amplitude modulation (PAM). However, the 
general architecture may be applied to other modulation techniques such as quadrature 
amplitude modulation (QAM), carrieriess amplitude and phase modulation (CAP), etc. 
The technique utilizing the identification algorithm is illustrated in Figure 12, 

10 depicting a signal flow of the joint co-channel identification and symbol detection 
architecture based on a batch identification algorithm. In Figure 12, y^Jn) is the 
aggregate disturber, s{n-5) is an estimate of the PAM symbols sent at time n-5, 5 is the 
delay introduced by the PAM receiver, a nd h is a vector that contains the co-chan nel 
model parameters. If L disturbers are present in the aggregate disturbance >^^^./nj, 

15 s(n - 5) is an L-dimensional vector. 

For explanation purposes, assunie that aniniU a^ of the co -£hanne 

impulse response, jip ^ prov ided. Then, using this initial model, we design a PAM 
receiver to detecuhe symbolsj^^ disturbance. When several disturbers are 

present in the mixture signal y^Jn), then the PAM receiver may be designed as a joint 

20 receiver, a successive receiver, or a parallel one. In any case, the PAM receiver can be 
selected from the well known variety of standard PAM receivers such as linear equalizers 
followed by a decision device, a decision feedback equalizer, a Viterbi algorithm, etc. 
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The selection of the receiver structure depends on the signal to noise ratio of the received 
signal, the amount of computational resources available, etc. 

The output of the PAM receiver is used as the estimated input in a batch 
identification algorithm. Similarly, the output corresponds to the aggregate disturbance 
5 appropriately delayed by 5. The batch identification may be a system identification 
algorithm as is well known to those in the art. A batch identification algorithm is 
illustrated in Figure 13. For purposes of explanation of this embodiment, we assume that 
a particular model structure has been previously selected. All the possible model 
structures may be grouped into two large categories: finitejm{ni^^ 

10 models andjnfinite impulse response (IIR). models. The selection of one structure or 
another may depend on the characteristics of the co-channels to be identified. The first 
stage of the algorithm is to collect N points of the input and output signals. Then, 
according to the model structure previously selected, the regression matrices are formed 
using the recorded inputs and outputs. Finally, the model parameters are obtained by 

15 solving a least squares problem. A number of computationally efficient least squares 
algorithms can be used to solve the problem like square root algorithms, QR 
factorizations, etc. 

If we let h, be the co-channel impulse respo nse obtained usin g the batch 
algorithm, then once the batch identification algorithm is completed, the switch in Figure 
20 p is switched from the initial co-channel position to the.ID'd co-channel position. The 
new^^demified^ used to re-design the PAM receiver and the procedure is re- 

iterated K times until convergence. The co-channel obtained after the ^-th iteration is h,. 
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Note that if several disturbers are present, this strategy may be combined with a 
successive interference cancellation algorithm to successively obtained models for the co- 
channels of the different disturbers. 

A second alternative for implementing block 614 in Figure 6 is to use an adaptive 
5 channel tracking technique such as the one shown in Figure 14. In this case, we 

implement a PAM receiver subsystem, 1410, and a parameter adaptation algorithm 1430. 
The particular details of these two blocks will be explained in a later section. 

Parameter Adaptation 

10 An initial assumption for initial co-channel identification relied upon an 

observation time short enough to neglect instantaneous phase errors of the timing signal 
due to frequency drift or other random effects. Thus, the initial identification results 
relied upon the baud rate being known and fixed. In some applications, it is important to 
track the changes of the co-channel impulse response during an extended period of time. 

15 In those cases, baud rate drift and random effects will affect the performance of the 
particular application over long time periods. 

According to ITU-T recommendation G.810 "Definitions and terminology for 
synchronization networks*', 08/96, the instantaneous phase of a timing signal is represented as 



20 where <l>o represents the initial phase offset, 5{t) is related to frequency deviation and drift, 
and (pit) corresponds to the random phase deviation including both clock jitter and wander. 
For the second term in the above equation, the combined effect of frequency deviation and 
drift could be on the order of 60ppm (parts-per-million). Variations caused by 5(0 and (p(t) 



cD(0 = Oo+5(0 + <p(0 
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pose a serious problem to several applications that require a long time observation period. 
Thus, the timing issue must be addressed in order for these applications to work effectively. 
This is accomplished in block 616 of Figure 6. 

What the present invention presents is a data-aided adaptive channel tracking 
technique to address this issue. In particular, we will describe the procedure illustrated in 
Figure 14. Block 1410 is a PAM receiver. As it has been previously indicated for Figure 12, 
the PAM receiver can be selected from the well known variety of standard PAM receivers 
such as linear equalizers followed by a decision device, a decision feedback equalizer, a 
Viterbi algorithm, etc. The purpose of the PAM receiver is to take as an input the aggregate 
disturbance and produce an estimated PAM symbol sequence. We can use those estimates in 
a data-aided adaptive algorithm 1430 to track the time-varying co-channels. A possible 
implementation for the adaptation algorithm is a recursive least squares (RLS) algorithm. 
Another implementation is a least-mean squares (LMS) algorithm. It will be appreciated by a 
person of ordinary skills in the art that a rich variety of adaptation algorithms can be 
implemented in block 1430. The process is illustrated in Figure 12, where y^i,,{n) is the 
aggregated disturbance signal corrupted by co-channel inter symbol interference (ISI), DMT 
residue and additive color Gaussian noise, s(n) is the estimated PAM symbol sequence, and 

h is the updated co-channel impulse response. 

The tracking procedure can be summarized as follows: 

(1) Initialize the co-channel impulse response with the one identified in the initial co- 
channel identification 608; 

(2) Run the PAM receiver and output the estimated PAM symbols s{n) ; 

(3) Run the parameter adaptation algorithm to obtain the estimated channel h ; 
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(4) Provide the estimated h to the PAM receiver for use in the next segment of data. 

Application of the present invention in one embodiment in an ADSL svstem 

Figure 15 illustrates an example of how the various parts of the identification 
5 process may prove useful in an interference compensation device incorporated in an 
ADSL modem. A typical sequence of events is shown in Figure 15 starting with initial 
power being applied to the modem at 15 IQ. Next, the modem will enter a training time 
1520, in which tasks such as time equalization (TEQ) training 1522 and frequency 
equalization (FEQ) training 1524 may occur. Please note that TEQ and FEQ training are 
10 standard within an ADSL modem and are shown here to help the reader appreciate and 
understand where the present invention fits within the overall picture for one 
embodiment. 

After the training time 1520 is completed, the next step is identification 1530 of 
possible crosstalk sources. Within the identification 1530, may be tasks such as detection 

15 of service types present 1532, baud-rate estimation 1534, setup of co-channel estimation 
1536 and initial co-channel estimation 1538. After identification 1530 of possible 
crosstalk sources has been completed, the next step is system design 1540. The system 
design 1540, may include, for example, such tasks as compensator design 1546,~and a 
final co-channel estimation 1548. For an example of compensator design see co-pending 

20 patent application Serial No. titled "Method and Apparatus for 

Mitigation of Disturbers in Communication Systems" assigned to the assignee herein and 
filed on even date herewith. 
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When system design 1540 of compensation is completed and deployed, the DSL 
modem may be used for customer-initiated communications at transmission time 1550. 

Transmission time 1550, also sometimes referred to as showtime, may include, 
compensation deployment 1554 and parameter adaptation 1556. For an example of 

5 compensation deployment see co-pending patent application Serial No. 

titled "Method and Apparatus for Mitigation of Disturbers in Communication Systems" 
assigned to the assignee herein and filed on even date herewith. 

Part of the modem parameter adaptation 1556 may be accomplished via the 
channel tracking procedures described in the current invention. Finally, if the modem is 
10 turned off we have an end 1560 the operation. The present invention discloses 

techniques that are applicable to various blocks of Figure 15, with major emphasis on 
Identification 1530, Final Co-Channel Estimation 1548, and Parameter Adaptation 1556. 

Note that while the operations of Figure 15 are illustrated as sequential steps, this 
is not the only embodiment possible. For example, identification 1530 may involve 
15 several processes such as detection of service types present 1532 and, for example, baud- 
rate determination occurring concurrently. Likewise, while at transmission time 1550, a 
compensation system may be deployed as completed at system design 1540, there is 
nothing to preclude rnore identification 1530 of crosstalk sources during transmission 
time 1550. That is, for example, the identification 1530 may be a batch mode 
20 identification or periodically invoked, or even continuous in nature. 

Thus, a method and apparatus for identification of interference, such as crosstalk 
sources, and their coupling channels are disclosed. Although the present invention has 
been described with reference to specific exemplary embodiments, it will be evident that 
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various modifications and changes may be made to these embodiments without departing 
from the broader spirit and scope of the invention as set forth in the claims. Accordingly, 
the specification and drawings are to be regarded in an illustrative rather than a restrictive 
sense. 
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